GLOBAL ATTRACTORS AND DETERMINING MODES FOR THE 
3D NAVIER-STOKES-VOIGHT EQUATIONS 
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Abstract. We investigate the long-term dynamics of the three-dimensional Navier- 
Stokes-Voight model of viscoelastic incompressible fluid. Specifically, we derive upper 
bounds for the number of determining modes for the 3D Navier-Stokes-Voight equations 
and for the dimension of a global attractor of a semigroup generated by these equations. 
Viewed from the numerical analysis point of view we consider the Navier-Stokes-Voight 
model as a non-viscous (inviscid) regularization of the three-dimensional Navier-Stokes 
equations. Furthermore, we also show that the weak solutions of the Navier- Stokes- 
Voight equations converge, in the appropriate norm, to the weak solutions of the inviscid 
simplified Bardina model, as the viscosity coefficient v — > 0. 
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1. Introduction 

We consider the three-dimensional Navier-Stokes-Voight (NSV) system of equations 

v t - vAv - a 2 Av t + (v ■ V> + Vp = f(x), x G Q, t G R + , (1.1) 
div v = 0, xeQ,te R + ; v{x, t) = 0, x G dQ,t G K+, (1.2) 
v(x, 0) = vo(x), xefl, (1.3) 
where Q C M. 3 is a bounded domain with sufficiently smooth boundary dQ, v = v(x,t) is 
the velocity vector field, p is the pressure, v > is the kinematic viscosity, a is a length 
scale parameter characterizing the elasticity of the fluid, and / is a given force field. 

The system (]l.ip ~ (jl.2|) models the dynamics of a Kelvin- Voight viscoelastic incompress- 
ible fluid and was introduced by A. P. Oskolkov in [38] as a model of motion of linear, 
viscoelastic fluids. 

The viscous simplified Bardina model was introduced and studied in [34] (see also [1]) 
as a simplified version of the Bardina sub-grid scale model of turbulence [3]. In [3] the 
viscous and inviscid simplified Bardina model were shown to be globally well-posed. It is 
interesting to observe that the inviscid simplified Bardina model coincides with the invis- 
cid version of the NSV equations fll.ip - fll.3l) . Viewed from the numerical analysis point 
of view the authors of [5] proposed the inviscid simplified Bardina model (or equivalently 
the inviscid NSV equations) as a non-viscous (inviscid) regularization of the 3D Euler 
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equations, subject to periodic boundary conditions. Motivated by this observation the 
system fll.ll) - fll.3p was also proposed in [5] as a regularization, for small values of a, of 
the 3D Navier-Stokes (NS) equations for the purpose of direct numerical simulations for 
both the periodic and the no-slip Dirichlet boundary conditions. 

In [38] it is shown that the initial boundary value problem fll.ip - fll.3l) has a unique 
weak solution. In [25] and [26] it is shown that the semigroup generated by the problem 
fll.ip - fll.3l) has a finite dimensional global attractor. 

In this paper we give an estimate of the fractal and Hausdorff dimensions of the global 
attractor of a dynamical system generated by the problem fll.ll) -( fTT3l) . which is an im- 
provement of the estimates done in [26J. Moreover, we derive estimates for the number of 
asymptotic determining modes of the solutions of the problem f ll.ip - fll.3l) . We also show 
that there exists a number m such that each trajectory v (t) on the global attractor of 
the dynamical system generated by this problem is uniquely determined by its projection 
P m v{t) onto the spanjwi, w m } of the first m eigenf unctions of the Stokes operator. 
This observation is related to the notion of continuous data assimilations as it has been 
presented in [29j,[36j and [37J. 

It is worth stressing that by adding the regularizing term (— aAv t ) to the NS equations 
the system f ll.ll) - fll.3p changes its parabolic character. In particular, the 3D system fll.ip - 
fll.3p is globally well-posed forward and backwards in time. The semigroup generated 
by the problem fll.ll) - fll.3p is only asymptotically compact. In this sense the system is 
similar to damped hyperbolic systems. We also remark that this type of inviscid regu- 
larization has been recently used for the two-dimensional surface quasi-geostropic model 
[28] . In particular, necessary and sufficient conditions for the formation of singularity 
were presented in terms of regularizing parameter. 



2. Preliminary 

In this paper we will be using the following standard notations in the mathematical theory 
of NS equations: 

• L P (Q), 1 < p < oo, and H S (Q) are the usual Lebesgue and Sobolev spaces, respec- 
tively. 

• For v = (vi,V2,vs), and u = (ui, 112,113) we denote by 



3 3 3 

2 

'i^'lll, 2 ^)- 

j=l j=l j,i=l 



[ U ,v) = ^(Vj,Uj)za (n ), INI 2 = IHli 2 (n)> W VV W 2 := W d i V 



We set 

V := {v G (C °°(fi)) 3 : V-v = 0}. 
H is the closure of the set V in (L 2 (fi)) 3 topology. 

P is the Helmholz-Leray orthogonal projection in (L 2 (Q)) 3 onto the space H, and 
h:=Pf. 
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• A := —PA is the Stokes operator subject to the no-slip homogeneous Dirichlet 
boundary condition with the domain (H 2 (Q)) 3 n V. The operator A is a self- 
adjoint positively definite operator in H, whose inverse A~ x is a compact operator 
from H into H. Thus it has an orthonormal system of eigenfunctions {wj} c *L l of 
A. 

• We denote by {Xj}'jL l , < Ai < A2 < • • • , the eigenvalues of the Stokes operator A 
corresponding to eigenfunctions {wj}° : L l , repeated according to their multiplicities. 

• V s := D(A S / 2 ), \\v\\ s := \\A s / 2 v\\,s G R. V : = Vi = (#o(^)) 3 n H is the Hilbert 
space with the norm ||i>||i = ||w||v = ||Vtt||, thanks to the Poincare inequality 
P - Clearly Vb = H. 

• For u,v,w G V we define the following bilinear form 

B(u, v) := P ((u ■ V)f ) and the trilinear form b(u, v, w) = (B(u, v),w). 

The bilinear form B(-,-) can be extended as a continuous operator B : V x V — > V, 
where V is the dual of V (see, e.g., [XT]). 

• For each u, v, w e V 

6(tt, f,f) = 0, and 6(tt, f , w) = — w, f ). (2.1) 

Next we formulate some well known inequalities and a Gronwall type lemma that we 
will be using in what follows. 



Young's inequality 

"i/(p-i) 



ab < -a p H — -^b g , for all a,b,s > 0, with q = p/Op - 1), 1 < p < 00. (2.2) 

p qe L/{p L > 



Poincare inequality 

HI < Ai 1/2 ||w||i, Vm G V, (2.3) 

where Ai is the first eigenvalue of the Stokes operator under the homogeneous Dirichlet 
boundary condition. 

Hereafter, C will denote a dimensionless scale invariant constant which might depend on 
the shape of the domain Q. 

Ladyzhenskaya inequalities ( [TT] . [31] . [33] ) 

\\u\\ L s < l7||m|| 1/2 ||Vm|| 1/2 , Vm G V, (2.4) 

1Mb < l7|| m || 1/4 ||m||i /4 , Vm g V. (2.5) 
Sobolev inequality (see, e.g., [TJ) 

||w|U« < CIMIi, Vm G 1/. (2.6) 
Gagliardo - Nirenberg inequalities (see, e.g., [2]. [IT]. [33]) 

||w|U6/(3-2e) < C||u|| 1-e ||u||f, 0<e<l, Vm G V. (2.7) 
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\\u\\ LP < CWuW^WuW 1 -^, pe[2,oo), Vne% (2.8) 
Agmon inequality (see, e.g., [TT] ) 

\\u\\ L oo m < CWuWl^WAuW 1 ' 2 , \/u e V 2 . (2.9) 
We will use also the following estimates of the trilinear form b(u, v, w) which follow from 
(JH - m> (see, e.g., [11]). 

\b(u,v,w)\ < L7||n|| 1/2 |H|j /2 |H|i|H|i, Vu,v,weV, (2.10) 
\b(u,v,u)\ < CIHI^IMI^IHl!, Vn,nGV, (2.11) 
\b(u,v,w)\ < ClwIlilMdlMI^IMlj 72 , Vu,v,weV, (2.12) 

\b(u, v , w)\ < CA 1 ^ 4 ||w||i||t'||i||iy||i, Vu,v,w &V. (2.13) 

Lemma 2.1. ( |23] . see also [15] j Let a(t) and b(t) be locally integrable functions on (0, oo) 
which satisfy for some T > the conditions 

^ ft+T ^ rt+T ^ rt+T 

liminf— / a(r)dr = 7, limsup— / a~(r)dT = T, liminf— / b + (r)dT = 0, 

*->oo T J t t^oo T J t t^oo T J t 

where 7 > 0, T < 00, a~ = max{— a, 0} and b + = max{o, 0}. // a non-negative, absolutely 
continuous function <p(t), satisfies 

<f)'(t) + a(t)<j>(t) < b(t), t G (0, 00), 

then (p(t) ^ as t — > 00. 



Definition 2.2. (see, e.g., [IS], US], [22]) A semigroup S(t) : V -> V,t > is called 
asymptotically compact, if for any sequence of positive numbers t n — > 00 and any bounded 
sequence {t>„} C V the sequence {5(t n )f n } is precompact in V. 

Theorem 2.3. (see, e.g., [19], [32], [13]^ Assume that a semigroup S(t) : V — > V, /or 
£ > £o > can &e decomposed into the form 

S(t) = Y(t) + Z(t), 

where Z(t) is a compact operator in V for each t > i > 0. Assume also that there is 
a continuous function k : [to, 00) x M + — > M + sitcA tnat /or every R > fc(t, i?) — > as 
t — > 00 and 

< fc(* 3 -R), for all t>t >0, and all \\v\\ v < R. 
Then S(t) : V —> V,t > is asymptotically compact. 

Next we state a result from [32] which will enable us to estimate the dimension of the 
global attractor for the system fll.ip - (ll.3l) . This result is typically useful in the context of 
nonlinear damped hyperbolic systems, when the damping term is not strong enough to 
control the instabilities rising from the perturbed nonlinearity. 
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Theorem 2.4. (see [13], [32],) Let S(t),t G R + ; fre a semigroup generated by the problem 

in i/ie phase space H and let M. a H is a compact invariant subset with respect to S(t). 
Let S(t) and $(•) be uniformly differentiable on A4 and let L(t,Vo) be a differential o/$ 
at the point S(t)vo,vo G M.. Suppose that L c (t,v ) := L(t, vq) + L*(t, vq), t>o G M. satisfies 
the inequality 

m 

(L c (t)u,u) < -h (t)\\u\\ 2 + J2K(t)h\\l, (2.14) 

k=l 

for some numbers < 0, (k = 1, ...,m) and some functions h , h Sk G Ly oc (R), h Sk (t) > 

0,/io(t)>0/orai/teR + . 

T/ien 

dim H {M) < dim f {M) < N, 

where N is such that 

m 

-h (T) + J2K(T)N s *<0, 

for some T > 0. i/ere /ij(T) := ^ hi(r)dr. 

3. Existence of Global Attractors 

Applying the Helmholtz - Leray projector P to the system (jl.ip - (jl.2D . we obtain the 
following equivalent functional differential equation 

Vt + uAv + a 2 Av t + B(v,v) = h, h = Pf, (3.1) 

v(0) = v . (3.2) 
The question of global existence and uniqueness of (13.ip - (l3.2p first was studied in [38] . 
where actually it was established that the problem fll.ip - fll.3B generates a continuous 
semigroup S(t) : V — ► V,t G M + . In [S] the authors proved also the global regularity for 
inviscid model of (I3.ip . i.e. when v = 0. 

In this section we show that the semigroup S(t) generated by the problem fll.ip - fll.3l) 
has an absorbing ball in V and an absorbing ball in V2. Then we show that S(t) : V — > V, 
for t G M + is an asymptotically compact semigroup, and deduce the existence of a global 
attractor in V. 

Let us note that the formal estimates we provide below can be justified rigorously by 
using a Galerkin approximation procedure and passing to the limit, by using the relevant 
Aubin's compactness theorem as for the NS equations ( see, for example, [11], [15], [H] or 

m)- 

Absorbing ball in V. Taking the inner product of (13. ip with v, and noting that due to 
( JXID (B(v,v),v) = 0, we get 

4 [Ht)\\ 2 + ^\\v(t)\\l] + 2v\\v(t)\\l < 2|NI-i|Nt)||i. (3.3) 
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It is easy to see by Poincare inequality (I2.3P that 

HK*)II! > ~ [AilMI 2 + > [IK*)II 2 + ^\\v(t)\\l] , 

where d := | min{^2 , Ai} = vd\. Hence (13 .3p implies 

j t [Mt)\\ 2 + a 2 \\v(t)\\l] + d [\\v(t)\\ 2 + a 2 \\v(t)\\l] < ±\\h\\\. 
By Gronwall's inequality we have 

\Wt)\\ 2j rO?\\v{t)\\\< 

1 



e -d (t-s) 



\\h\\ 2 

|2 i 2|| / \||2 ll'll-l 



\v (s) + a f (s 



11 7~ 



2 1 + 4 W1 " {E > 



Therefore, 

limsup [|k(i)|| 2 + « 2 |kWlli] < ^T 1 - (Ei) 

t^oo vdo 

The last inequality implies that the semigroup S{t) : V — > V, t G M + generated by the 
problem (ll.ll) - (ll.3p (or equivalently (I3.1I) - (I3.2I) ) has an absorbing ball 



B x :=<veV: |M|i < -== \\h\U \ . (3.4) 
[ \Jva l do ) 

Hence, the following uniform estimate is valid 

\\v(t)\\i < M u (3.5) 
where Mi = f° r t large enough ( t ^> 1) depending on the initial data. 

Asymptotic compactness. By using the Galerkin procedure it is not difficult to prove 
the following 

Proposition 3.1. Let sGl. If Wo E V s , g E L 2 ([0, T); V s -q) then the linear problem 

z t + a 2 Az t + vAz = g(t), z{0) = (3.6) 

has a unique weak solution which belongs to C([0, T); V s ) and the following inequality holds 

sup \\z(t)\\ s < C\\g\\ L 2^ T . Vs _ 2) , sGi 
te[o,T) 

Proposition 3.2. Let h e H , be time independent, then the semigroup S(t),t > is 
asymptotically compact semigroup in V . 



Proof. Let vq G V. First we observe that S(t) has the representation 

3{t)v = Y(t)v + Z(t)v , (3.7) 
where Y(t) is the semigroup, generated by the linear problem 

y t + vAy + a 2 Ay t = 0, y(0) = v , (3.8) 
and z{t) = Z(t)(vo) is the solution of the problem 

z t + vAz + a 2 Az t = h- B(v(t),v(t)), z{0) = 0, (3.9) 
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where v is the solution of (]l.ip ~ (jl.3|) (or equivalently (I3.1l) - (l3.2p ) with the initial data v . 
Taking the H inner product of (I3.8P with y we obtain 

j t [WW + o?\\y{t)\tf\ + d [\\y(t)f + a 2 \\y(t)\\l] < 0, 

where we recall that do = vd\ = ^|min{^2, Ai}. 
This inequality implies that 

\\y(t)\\ 2 + a 2 \\y(t)\\l < e"** [IN| 2 + a 2 |N| 2 ] , for all t > 0. (3.10) 
So the semigroup Y(t) : V — > V is exponentially contractive. 

Due to Holder's inequality and the Sobolev inequality (12.61) we have 

11^(^,^)11-1/2= SUp b(v,V,(f)) = 



4>&V,\\A 1 / A 

sup / P((v ■ V)f ) • (j)dx = sup / (v ■ V)v ■ P(f)dx = 
0eV,||AV 4 <£||=i Jn ^evjjAVi^i^i Jn 

sup / (v ■ V)v ■ (j)dx < C sup IMU 6 IMIi 

\A 1 / i 6\\=i Jn 6GV.\\A 1 / i 6\\=i 



L3. 



4>&v,\\A 1 / i (j>\\=i Jn 4>ev,\\A 1 / 4 4>\\=i 
Hence due to the Sobolev inequality ||0||l 3 < CH^^VII an d (12.61) we have 

||£M)||-i/2 < C sup WvWlWA 1 / 4 ^ < C\\v\\l, (3.11) 

^V,\\A l / 4 '4>\\=l 

and 

B(v,v) GL°°(M + ;^ 1/2 ). 
The function v{t) as a solution of the problem (I3.ip - (I3.2I) with t> G V belongs to 
L°°(R + ;l / ). Thus due to the inequality (13.111) and the Proposition 13.11 the solution 
of the problem (I3.9P belongs to C(R + ; V3/2), that is the operator Z{t) maps V into V3/2. 
Since the embedding V3/2 C V is a compact embedding, the operator Z(t) is a compact 
operator for each t > 0. Hence, the semigroup S*(t) satisfies the conditions of the Theorem 
12.31 an d is an asymptotically compact semigroup. □ 

Since each bounded dissipative and asymptotically compact semigroup possesses a com- 
pact global attractor (see, e.g., [2], [T5] , [31] . [13] ) we have: 

Theorem 3.3. If h & H then the semigroup S(t) : V — > /ias an absorbing ball B\ = 
{v G V : || f ||i < Mi} and a global attractor A± G V . The attractor Ai is compact, 
connected and invariant. 

Next we show that the global attractor A\ is a bounded subset of V2. 

Taking the inner product in Vi/ 2 of the equation (13.91) with z, and remembering that 
v(t) = y(t) + z{t) G Ai, we get 

|[lkWlli/2 + « 2 lkWll3/ 2 ]+2HkWll3/2 = 

2(h, z(t)) 1/2 - 2(B(v(t), v(t)), z(t)) 1/2 . (3.12) 
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The first term on the right-hand side has the estimate 

\2(h,z(t)) 1/2 \ < 2||/ i ||_ 1/2 |K*)|| 3/2 < P -\\z{t)\\l /2 + hh\\\ /2 . 
The second term, due to (13.111) . has the following estimate 
\2(B(v(t),v(t)),z(t)) l/2 \ < C\\(B(v(t),v(t))\\-i/2Mt)\\3/2 < 

+ ^\\B(v(t)xt))\\ 2 -i/2 < ^\\m\t /2 + %\\v 

Taking into account the last two inequalities in (13.121) we obtain 
d 



, t l\\zm\/2+C?\\z(t)\\l,2] +2d Q [\\zm\/2+^Ut)\\l/2] < 

c 

""" Ml n 



h\\ 2 - 1/2 ) • 



-(n\\l /2 < 1 ^ Y - i {Mt + \\hf_ x/2 ) = Lq. (3.13) 



V 

Integrating the last inequality we obtain the estimate 

C_ 

Since the attractor A\ is invariant, S(t)Ai = At, and due to (13 . 1 0[) the inequality 

Hi) - zm, = um, < c(\\ y mi)e- dot 

holds, we deduce that for each u G A\ there exists a sequence {z(tk)},tk — > oo, corre- 
sponding to ffe(O) G A\, such that 

u = lim z(t k ), D fc (0) G Ai. (3.14) 

fc^oo 

Thanks to (13.131) the sequence {z(tk)} is belonging to a ball in V3/2, whose radius L 
depends only on Mi and \\h\\. Hence, the sequence {z(tk)} is weakly compact in V 3 / 2 . 
Thus, by using (13 . 141) and the inequality IHI3/2 < liminf^^oo 1 1 fc ) 1 1 3 /a , we see that Ai 
is bounded in V3/2. 

Knowing that Ai is bounded in V3/2 we can use similar arguments to show that Ai is also 
bounded in V5/3 and in V 2 . 

V 2 absorbing ball. To show that the semigroup S(t) : V 2 — > V 2 has an absorbing ball in 
the phase space V 2 = D(A) we take H inner product of (13.11) with Av(t): 

^ [11^(^)11? + « 2 ||^(t)|| 2 ] + 2^||^(^)H 2 -h 2(^(^(t),^(t)),^(t)) = 2(/z,^(t)). (3.15) 

For the first term in the right hand side of (I3.15P we have 

\2(h,Av(t))\ < ^\\h\\ 2 + u\\Av(t)\\ 2 . (3.16) 
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By using the Agmon's inequality (12. 9p and Young's inequality (12. 2p with p = 4/3 we can 
estimate the last term in the left-hand side of (13. 15j) as follows 



2\{B{v,v),Av)\ < tf||u||£» ( n)||||v||i||||Ai;|| < C"IMI? /2 || \\Avf /2 < 



-AM\ 2 + -M\i 



Employing (13.161) and the last inequality, with e = 2z//3, we obtain from (13.151) 

j f [\\v(t)\\l + a 2 \\Av{t)\\ 2 ] + "\\Mt)\\ 2 < l\\h\\ 2 + ^Ht)\\t- (3-17) 
It follows from (13.171) that 

j t [Ht)\\l + * 2 \\Mt)\\ 2 } + do [Ht)\\l + a 2 ||^(t)|| 2 ] <l\\h\\ 2 + ^|| 

Let t be so that (13.51) holds for all t > t . Then integrating the last inequality over the 
interval (to,^) we get 

\\v(t)\\l + a 2 \\Av(t)\\ 2 < 

[\\v(t )\\l + « 2 ||^(t )|| 2 ] e- d ^-^ + ^ (1 - e- d °^) , (3.18) 
where R 2 := l\\h\\ 2 + §Mf. 

The last inequality implies existence of an absorbing ball 

B 2 :={veV 2 :\\Av\\<M 2 }, (3.19) 
where Mf = (ct2 ffi )do - That is, for all t » 1, we have \\Av(t)\\ < M 2 . 

Similarly, we can prove the following theorem 

Theorem 3.4. If h G V\, then the semigroup S(t) : V 2 — > V 2 has a global attractor 
A 2 C V 2 . The attractor A 2 is compact, connected and invariant. Moreover, A 2 is a 
bounded set in V3. 

Remark 3.5. Let us note that in case we assume in Theorem 13.31 that h G V±, instead of 
h G H, then the attractors A\ and A 2 coincide. 

4. Estimates for the Number of Determining Modes 

It is asserted, based on physical heuristic arguments, that the long-time behavior of tur- 
bulent flows is determined by a finite number degrees of freedom. This concept was 
formulated more rigorously for 2D NS equations by introducing the notion of determining 
modes in [17] . In [T7J it was shown that there exists a number m such that if the first m 
Fourier modes of two different solutions of the NS equations have the same asymptotic 
behavior, as t — > 00, then the remaining infinitely many number of modes have the same 
asymptotic behavior. 
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In [31] it was shown that the semigroup generated by the initial boundary value problem 
for the 2D NS equations with Dirichlet boundary condition has a global attractor which 
is compact, invariant and connected . It was also established in [31] that there exists a 
number m such that if projections of two different trajectories on the attractor on the m 
dimensional subspace of H, spanned on the first m eigenfunctions of the Stokes operator, 
coincide for each t 6 R, then these trajectories completely coincide for each t G R. 

The results obtained in [T7] and [31] were developed, generalized, and applied to var- 
ious infinite dimensional dissipative problems (see, e.g., [7].[8].[^.[T5].[T6].[T8].[20].[22]. 

and references therein). 



In this section we are going to give estimates for the number of determining modes 
(both asymptotic and for trajectories on the attractor) for 3D NSV equations. 

Asymptotic determining modes. Let us denote by P m the L 2 -orthogonal projection 
from H onto the m- dimensional subspace H m = span {w\, W2, w m }. We set Q m = 
I — P 

Let v and u be two solutions of NSV equations 

v t + vAv + a 2 Av t + B(v, v) = h(t), v(0) = v , (4.1) 
u t + vAv + a 2 Au t + B(u, u) = g(t), v(0) = v . (4.2) 

Definition 4.1. A set of modes {w\, ■ ■ ■ ,w m } is called asymptotically determining (see 

mm) if 

lim \\v(t) -u(t)||i = 

whenever 

lim \\h(t) - g{t)\\-i = and lim \\P m {v(t) - u(*))||i = 0. 

t— >oo t^oo 

Theorem 4.2. Assume that the following conditions are satisfied 

|_i < h < 00, Vt e R. (4.3) 



lim \\h(t) - 0(t)||_i = and lim \\P m (v(t) - u(t))\\ = 0. (4.4) 

t— KX) t^OO 

Then the first m eigenfunctions of the Stokes operator are asymptotically determining for 
the NSV equations with homogeneous Dirichlet boundary conditions, provided m is large 
enough such that 

A ™ +1 > c ^ (4 - 5) 

Proof. It is clear that the function w = v — u satisfies 

w t + vAw + a 2 Aw t + B(v, w) + B(w, v) - B(w, w) = 9(t), v(0) = v , (4.6) 
where 6(t) = h(t) -</(*). 
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It is clear from the proof of (Ei) that 

h 

limsup ||v(t)||i < 7=. (4.7) 

Multiplying (14.61) by q(t) = Q m w(t) in H we obtain 

j t [\\q\\ 2 + a 2 \\q\\l]+2v\\ q \\l + 2b( q ,v,q) = 

2(9, q) - 2b(v, p, q) - 2b(p, p, q) + 2b(q, p, q), (4.8) 

where p = P m w. 

Before estimating the terms of (14.81) we observe that for each <p G V we have 

||Q m 0||i > A m+1 ||Q m 0|| and ||P m 0||i < A m ||P m 0||. (4.9) 
Due to the inequality ( 12. lip the term b(q, v, q) has the following estimate: 

2\b(q,v,q)\ < CUnqf/'MU < C^v^. (4.10) 

A m +i 

The first term in the right-hand side of ( 14.81) has the estimate 

m?)l<^llCi + ~IMIi- (4.ii) 

Employing the inequalities (I2.12p and (14.91) we estimate the second term in the right- 
hand side of ( 14. 8 p as follows 

2\b(v,p,q)\ < ||^||i|b||i|kH 1/2 |k||i /2 < CX m X- l _H\\ph (\\q\\j + \\v\\l) . (4.12) 

Other terms in the right-hand side of (14.81) can be estimated in a similar way to (I4.12p . 
Using estimates (I4.10p - (l4.12p and the estimates of other terms in the right-hand side of 
(@~3D we obtain 

j t [H? + « 2 |Mli] + ^IMIi + \\q\\l (u--^ r \\v\\ 1 \< b(t), (4.13) 



A m+1 



2|| l|21 i V II ||2 i II ||2 / ' 

+ " IklliJ + + IkL v ~ ~i 

where b(t) is satisfying the corresponding condition of Lemma 1. 

Let us choose t\ > so large that ||v(£)||i < M\, for all t > t\ and m so that n(m) : = 
A m +i — (^t 1 ) 4 > 0. Then it follows from the last inequality the following relation 

J t [h\\ 2 + ^\\q\\i]+^\\q\\i<b(t), foraUt>t lf 

or 

j f [\\l\\ 2 + « 2 |klli] + d m [IMI 2 + a 2 ||g|| 2 ] < bit), for all t > h, (4.14) 
where d m = ~ min{^ , A m+i }. 

Thus, due to Lemma 1 the statement of the theorem follows. □ 
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Remark 4.3. Let us observe that the number m, for which A m+ i > holds, is an 
upper bound for the minimal number of asymptotically determining modes for weak 
solutions (i.e., solutions belonging to L°°(R + ;H) D L loc (R + ; V)) of the initial boundary 
value problem for the 3D Navier Stokes equations. 
In fact, for weak solutions of NS equations instead of (14.131) we have 

j t h\\ 2 + ^i (^Zi-cWviu) \\ q \\*<b(t), 

and instead of (14. 7\i we have for weak solutions of NS equations (see, e.g., [TT] . [TP] , [20] 
and [131) 



1 r +1 

limsup- / \\v(r)\\ldT < 

t-^oo J- Jt 



Tz/ 3 A? i/2 Ai 
Hence 

1 f t+T h h 

limsup— / ||u(r)||idT< 



Thus, the function a(t) := X^+i ^^m+i — ^IMIi J satisfies conditions of Lemma |2~TT1 
provided T is large enough and 

h 4 

A m +1 > C o , n • 

Different estimates of asymptotic determining modes for weak solutions of 3D NS equa- 
tions are obtained in [12] (see also [H] , [15] and references therein) . The estimate obtained 
in [12] involves generalization of the so called mean rate dissipation of energy, per mass 
and time, i.e. it involves 



l y* 

v limsup- / sup || Vf(x, r) || 2 dr. 



For other related results concerning estimates of the number of asymptotic determining 
degrees of freedom for weak solutions of the 3D NS equations see, e.g., [10], [20] and 
references therein. 

Determining modes on the attractor. Next we give an estimate of determining modes 
for trajectories on the attractor. 

Definition 4.4. A set of modes {w\, ■ • ■ , w m } is called determining on the attractor (in 
the sense of [31]) if for each two trajectories v (t) and u(t) on the attractor A\ the equality 

\\P m (v(t) -u(t))\\i = 0, foralUeR 

implies 

v(t)=u(t), VteR. 

Let v and u be arbitrary two trajectories in the attractor A\ of (3.1). Then w = v — u 
satisfies 

w t + a 2 Aw t + vAw + B(w, v) + B(u, w) = 0. (4.15) 
Taking the inner product of H4. 1 51) with q = Q m w we get 

J t [Ikf + IMI?] + MqWI = - 2b (™, V, q) - 2b(u, w, q). (4.16) 
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Assume that P m w(t) = 0, for all t £ R, then Q m u> = q satisfies 

j t [\\q\\ 2 + \\q\\l] + 2u\\q\\l = 2b(q,v,q). (4.17) 

Due to (12. lip we have 

|26(g,u,g)| < C||g||} /2 ||g||J /2 |H|i 
Noting that on the attractor A\ we have ||w||i < M\, we employ the last inequality, and 
inequality (14.91) to obtain from (I4.16P 

j t [Ikll 2 + « 2 |kll?] + u\\q\\l + |kir /2 ||g||f (uXHU - CM,) < 0. (4.18) 
Let us choose m, large enough, so that X m +i > (^r^) 4 - Then (I4.18P implies 

|[ll# + « 2 W+^[NI 2 + a 2 ||g|| 2 ]<0, 
where / m = |min{A m+1 , ^}. 

Finally, we integrate the last inequality and get 

||g(t)|| 2 + a 2 ||g(t)|| 2 < exp[-/ m (t - s)) [\\q(s)\\ 2 + a 2 \\q(s)\\l] . (4.19) 
Passing to the limit as s — > — oo we obtain 

||g(t)|| 2 + a 2 \\q{t)\\\ = 0, for all t E R. 
Thus, the following theorem is true. 

Theorem 4.5. Let v andu be two solutions of the problem U.1\) - I[T~3\) from the attractor 
Ai. Assume that P m (u(t)) = P m (v(t)),Vt £ R ; where m is so that 

A ™ +1 a C 5W (4 ' 20) 

Then v{t) =u(t), for allt £ R. 

5. Estimates of Dimensions of the Global Attractor 

In this section we show the differentiability of the semigroup with respect to the initial 
data. This is to prepare for implementing Theorem 12 .41 in order to estimate the dimension 
of the global attractor. 

Theorem 5.1. Let uq and vq be two elements of V. Then there is a constant K = 
K(\\uq\\i, ||vo||i) such that 

\\S(t)v -S(t)u -A(t)(v -u )\\i < K\\v -u \\l (5.1) 

where the linear operator A(t) : V —> V, for t > is the solution operator of the problem 

& + c?Ai t + Ai + B(i, v) + B(v, - 5(£, £) = 0, £(0) = v - u , (5.2) 

and v(t) = S(t)v . That is, for every t > 0, the map S(t)v , as a map S(t) : V — > 
V is Frechet differentiable with respect to the initial data, and its Frechet derivative 
D Vo (S(t)v )w = A(t)w . 
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Proof. It is easy to see that the function T)(t) := v(t) — u{t) — £(£) = S(t)(v — u ) — £(t) 
satisfies 

T] t + a 2 Ai] t + uAr) + B(rj, v) + B(v, rj) — B(w, w) = 0, 
where w = v — u. Taking the inner product of the last equation with 77 we obtain 
d 



(lf liriu +a 2 \\ V \\l]+2u\\ V \\l = -2b( V ,v, V )-2b(w,w, V ). (5.3) 

By using inequalities (13.51) and (I2.5P and Young's inequality we can estimate the terms in 
the right-hand side of (I5.3P as follows. 
By ( EH]) we have 

\2b( V ,v,r,) I < C\\vU V \\ 1/2 \\v\\l /2 < CMMvimi 3 / 2 < ^(IM| 2 + 3|M| 2 ). 



By Q23DD 



126(^,^,^)1 = 126(^,^,^)1 < CA^IMIilhH? < HMIi + T^rlMli- 



4i/Ai 



Hence, we obtain from (15. 3p 



j t dmi 2 + « 2 imi 2 ] < ^p(Ni 2 + mil) + 4^\M\i. (5-4) 

The function w(t) = v(t) — u(t) = S(t)v — S(t)u satisfies 

w t + a 2 Aw t + uAw + B(w, v) + B(v, w) — B(w, w) = 0, w(0) = v — u := w . 
Taking the inner product of the last equation with w, and using (I2.13P and (E2) we obtain 

j t [\\w\\ 2 + a 2 \\w\\\] + 2v\\w\\\ = 2b(w,v,w) < 2CA} /4 HilMli < 

[|M| 2 + a 2 |M|i] , 

where h = 2CAi / V 4 [\\v{0)\\ 2 + a 2 \\v{0)\\ 2 + (l/vd ) \\h\\ 2 _ J 1/2 . Integrating last inequal- 
ity we get 

lkWlli<(l + T^a)IKO)lliexp(«iO. (5.5) 
Aicr 

It follows from (j5~^|) and (1531) that 

I [hf + « 2 |l^l| 2 ] < M [H? + « 2 |M| 2 ] + A 2 \\w(0)\\i exp(2M). 
Integrating and using Gronwall's inequality : 

11^)11? <A(t)\\w(0) ||* 
where := exp(2/ti + Ai)t. So we have 
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Thus the differentiability of S(t) with respect to the initial data follows. 
We rewrite (13 .ip in the following form 

v t = ~v + ^G- 2 v - G- 1 B{G- 1 v, G^v) + G^h, (5.6) 
or a 2 

where G 2 = I + a 2 A, and v = Gv. 

The equation of linear variations corresponding to (15. 6p has the form 

w t = L(t)w, (5.7) 

where 

L(t)w := + ^G~ 2 w - G^BfG^w, G~ l v) - G- l B(G- l v, G~ l w). 

a 2 a 2 

Now we consider the quadratic form 

(L(t)w,w) = -- ^|M| 2 + ^WG^wf - biG^w, G-H.G^w). 
a 2 a 2 

By using inequality (12.111) and the inequality < we get 

\b(G~ 1 w, G~ l v, G^w)] < ^IG^HI^IMI^IHI. 

or' 2 

Employing Young's inequality with p = 4/3, e = 2z// (3a 2 ), and the fact that on the global 
attractor Ai the estimate ||u|| < (Ai + a 2 Y^ 2 Mi holds, we obtain 

IHG-V G-\ G^w)\ < ^-\\w\\ 2 + C(Xl+ ff Mt \\G-'w\\ 2 - 

la 2 i/^a 4 

Due to the last inequality the quadratic form (L(t)w,w) has the following estimate 

(L( t ) W , W )<-- M 2 + (- + v J ^ \\G-^\\ 2 (5.8) 

Thus, we can use Theorem 12.41 to get the desired estimate for the fractal dimension of the 
attractor Ai 



dMl) < c^pi + 2 < c (W)w, + 2 . (5 . 9) 



We recall that M\ = „J^ \\h\\-ii di — |min{a 2 , Ai}. Let us note that in our situation 



h (t) = —,s = 0, Sl = -1, /i ai t = — + 1 1 3 / 

2a 2 a 2 p a a^ 

and h Sk {t) = 0,k>2. 



□ 
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6. The Inviscid Limit 

Here we show that when v — > the weak solution of the initial boundary value problem 
for the NSV system, i.e. of the problem (ll.ll) - (ll.3p . is tending to the weak solution of the 
initial boundary value problem for the inviscid simplified Bardina model 

u t - a 2 Au t + (u ■ V)u + Vp = f, xGfi,t>0, , . 

V-u = 0, xeQ,t>0, ^ ' 

u(x, 0) = v (x), x G O; u(x, t) = 0, x G dQ, t > 0. (6.2) 

The problem of existence and uniqueness of solutions of the initial boundary value 
problem, with periodic boundary conditions, for the 3D viscous and inviscid simplified 
Bardina models is studied in [5]. In particular, it is shown in [5] that the problem ( 16. 1D - 
(16. 2p has a unique solution u G C-^R; V), for initial value u G V. 

Applying to (16. ip the Helmholtz-Leray operator P we obtain the equivalent functional 
differential equation 

u t + a 2 Au t + B(u,u) = h, (6.3) 

w(0) = v Q . (6.4) 

Let v(t) be the solution of (16. ip with initial v(0) = v G V. Denote by w = v — u. Then 
w satisfies the relation 

w t + a 2 Aw t + B(w, v) + B(u, w) = —uAv, (6.5) 

iw(0) = 0, (6.6) 

which holds in the space V. Taking the action of (16.51) on w, which belongs to V, and 
using a Lemma of Lions-Magenes concerning the derivative of functions with values in 
Banach space (cf. Lemma 1.2 Chap. III-p. 1 69- [44J ) , we obtain 

j- [\\w\\ 2 + a 2 \\w\\ 2 ] = -2v{Vv, Vw) - 2b(w, v, w). (6.7) 



For the first term in the right-hand side we have 

\2u(Vv, Vw)\ < v 2 \\v\\l + \\w\\l 
The second term we estimate by using the inequality (I2.13P 

\2b{w,v,w)\ < CXl^WvWxWwWl 
Utilizing last two inequalities in (16.71) we get 
d 



j t [\\w\\ 2 + a 2 \\w\\ 2 ] < v 2 \\v\\ 2 + (l + C\\'*\\v\\i) \\w\\l < 

+ a~ 2 (l + 2L7A} /4 ||i;|| 1 ) [||^|| 2 + a 2 |M|l] • 



u 2 \\vf 2 -- 2 1 1 -V'\ J Ml, ,n ^ in, ,n-' 

Integrating the last inequality and using the standard Gronwall's lemma we get the esti- 
mate 

].rll)\\ 2 + a 2 \\w(t)\\l < I |Kr)||?rfTexp|4 + ?^ /" ||Wr)||,./r I . IC,f 



a 2 



o 
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Next we show that on each finite interval [0, T] we can estimate ||w||i by a constant 
depending only on ||vo||, ||^o||i an d the parameter a. Indeed, (13.31) implies that 

j t [\Ht)f + a 2 \\v{t)\\\] <or*\\h\\\ + o?\\v{t)\\l 

Integrating the last inequality over (0, t) with respect to time variable we obtain 

\\v(t)\\ 2 + ® 2 Mt)\\l < lko|| 2 + « 2 ||^o||? + to- 2 ||/i|| 2 1+ f [\\v(r)\\ 2 + a 2 \\v(r)\\l}dT. 

Jo 

By using the Gronwall inequality we get 

IK*) Hi < D T e T , for alU G [0,T]. 
Here D T := a~ 2 [\\v \\ 2 + a 2 \\v \\ 2 + Ta- 2 \\h\\ 2 _ J . 

Hence (16.81) implies 

\\w(t)\\ 2 + a 2 \\w(t)\\ 2 < u 2 TD T e T exp (V 2 T + 2Ca- 2 \ 1 ' i TD l J 2 e T / 2 ) . (6.9) 

Remark 6.1. The problem of convergence of solutions of the NSV equations to solutions 
of NS equations as a — > was studied in [39J. It is shown in [5H] that strong solutions 
of the NSV equations converge to strong solutions of the NS equations as a — > 0, under 
specified smallness conditions on the initial data of the problem. 

Remark 6.2. The results obtained in this paper are valid also for the solutions of the 
initial boundary value problem for the 3D NSV equations with periodic boundary condi- 
tions. 

Finally we would like to notice that the results reported here can be extended to other 
similar equations, a subject of future work. For instance, for the 3D equations of motion 
of Kelvin- Voight fluids of order L > 1: 

L 

v t + {v- V)v - fx Av t - niAv - ^2 & Au i + = /, 

i=i 

dtui + aim - v = 0, l = l,...,L. 

where //o, Pi, ati > 0,1 = 1,...,L. Also for the generalized Benajamin-Bona-Mahony 
(GBBM) equation: 

u t -a 2 Au t + vAu + V ■ F(u) = h, (6.10) 
where a smooth vector field F(u) satisfies the growth condition 

\F(u)\ <L7(1 + M 2 ). 

The problem of existence of a finite dimensional global attractor and estimates for the 
number of determining modes on the global attractor of Kelvin- Voight fluids of order 
L > 1 is established in [27\. In [35] the existence of a finite dimensional global attractor 
is established for ID GBBM equation under periodic boundary conditions. Existence of 
a finite dimensional global attractor for 3D GBBM under periodic boundary conditions is 
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proved in j6]. In |42j it was shown the existence of the global attractor for GBBM equation 
in H 1 (R 3 ). Moreover, the existence of a global attractor for a similar two-dimensional 
model describing the motion of a second-grade fluid is established in |35j . 
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